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Abstract 

In a recent paper [Phys. Rev. B57, 3375 (1998)], we examined in detail 
the nonlinear (electrical) dc response of a random resistor cum tunneling bond 
network [RRTN, introduced by us elsewhere to explain nonlinear response of 
metal-insulator type mixtures). In this work which is a sequel to that paper, 
we consider the ac response of the i?i?TA^-based correlated RC {CRC) model. 
Numerical solutions of the Kirchoff 's laws for the CRC model give a power-law 
exponent (= 0.7 near p = p^) of the modulus of the complex ac conductance at 
moderately low frequencies, in conformity with experiments on various types 
of disordered systems. But, at very low frequencies, it gives a simple quadratic 
or linear dependence on the frequency depending upon whether the system is 
percolating or not. We do also discuss the effective medium approximation 
{EM A) of our CRC and the traditional random RC network model, and 
discuss their comparative successes and shortcomings. 
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1 Introduction 



Electrical response of a composite material made up of conducting and insulating 
components to an external field is still of much interest. To find the electrical 
conduction properties of random mixtures one usually models a network with the 
notions of percolative processes |^ as a basic framework. For simplicity one may 
consider a lattice-based percolation and assume that the occupied bonds (with a 
probability p) are metallic or good conductors and the rest of the unoccupied bonds 
(with a probability 1 — p) are insulators or very poor conductors. So it turns out 
that one is dealing with a random electrical network where each bond of the regular 
lattice constitutes an appropriate circuit element. To study the linear direct current 
(dc) response of a random binary mixture, one considers it to be a random (linear) 
resistor network {RRN) where the conductance (7 of a circuit element assumes the 
value Qi with probability p, or (72 (= for insulators) with the complementary 
probability 1—p. When one attempts to study the alternating current (ac) response 
of such a system at a frequency / = uj/{2tt), the elementary conductances gi and 
g2 do, in general, become complex quantities, termed as complex 'admittances' 
(inverse of impedances), because of the presence of inductances and/ or capacitances 
in some parts of the circuit. Electrical conductivities (both dc and ac) of various 
composite systems, have been the subject of numerous studies 0, |], |]. On the 
theoretical side, analytical calculations involving related Maxwell's equations with 
appropriate boundary conditions, effective medium approximation {EM A), series 
expansion method, Monte Carlo simulation, model random resistor network {RRN), 
etc. are employed to study various interesting behaviors. 

The key experimental facts to be contended with by a model in the case of non- 
linear dc response have been discussed in our recent works [^,0]. To reproduce the 
very low percolation threshold, a 'low'-temperature resistance minimum, and the 
nonlinear electric (dc) response of metal-insulator composites, we had introduced in 
our original work a new semi-classical bond percolation model. The distinguish- 
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ing feature in this model (compared to earlier RRN, etc.) is the possibility of a 
local tunneling/ dielectric breakdown (beyond a microscopic threshold voltage, Vg) 
between two metallic bonds separated by a nearest neighbor lattice distance. We 
call these specially positioned insulating bonds, tunneling or 't-bonds'. All the other 
insulating bonds separated farther away remain insulating at any voltage however 
large. Clearly, since the positions of these t-bonds are totally correlated to the posi- 
tions of the metallic bonds in a particular (random) configuration, this model may be 
considered as a correlated bond percolation model. The percolative phase transition 
aspects [H and the dielectric breakdown aspects for this model have already been 
discussed elsewhere. Based upon this correlated bond percolation, the circuit made 
of the random network of resistors and t-bonds has been called a RRTN (random 
resistor cum tunneling bond network). 

In the present work which is a sequel to our work on the nonlinear dc response of 
the RRTN (referred to as I from now on), we are concerned with the modeling 
of the generic ac response of a whole variety of experimental systems within the 
RRTN framework (reported briefly first in the ref.]^). Experiments on the complex 
ac conductance G^u) in various composite systems, dispersed metals etc., as well as 
many disordered/ amorphous systems 0, ^, ^ ||, |10|, |1T], |T^, report a non-integer 
power-law behavior of the modulus, |G'(co')| = [(-Re G)^ + (Jm G)^]^/^. At a fixed 
voltage and at a moderately low-cj, [|G(ct;)| — G(0)] oc cu", where a = 0.7 (the case 
of very low u's is discussed later). Now, at a low voltage V (see I), [G'(V) — G{V = 
0)] oc {V~VgY Q, where Vg is the macroscopic threshold voltage for dc nonlinearity 
and the exponent 6 is dimensionality-dependent. There has been a recent claim 
for 3D carbon- wax composites that 6 = a = 1.4. We shall comment on this later. 

In the traditional RC network models used to study the ac response of compos- 
ites, the conducting bonds are pure (real) resistors and all the insulating bonds may 
behave as capacitors in the presence of an ac-field. A fairly complete review and 
references (from the percolative aspect) on the linear response may be found in the 



ref. 0]. Another more recent review in this regard is due to Shalaev [|13 
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If one applies an ac electric field across our RRTN model, a t-bond between 
two nearby metallic bonds are expected to behave as a capacitor. Note that perfect 
capacitors at all the t-bonds correspond to a situation where all the t-bonds have 
zero dc conductance at low voltages, and hence the RRTN is in its lower linear dc 
(or, ohmic) regime. Similarly, /ea% capacitors with a very low constant conductance 
for all the t-bonds, implies that the RRTN is in its upper linear dc regime. To obtain 
the nonlinear ac response in the truly nonlinear dc regime, one has to let the t-bonds 
be active or passive according to the voltage differences across them. In this work 
we study the nonlinear ac response in either the upper or the lower linear dc regimes 
only. Further, we make a simplifying assumption that all the capacitances across 
insulators farther than the nearest neighbor distance are zero. Based on the RRTN 
model, this model for studying ac response may thus be called a correlated EC 
(CRC) model. We believe that the simplicity of our model is physically appealing 
and realistic enough. Since the capacitors placed in this way never percolates by 
themselves, the |G'(c(j)| is expected to connect nonlinearly between its own lower and 
upper ac saturation regimes of c<j ^ 0^ and — oo respectively. 

To make the above discussions more concrete we show in Fig. 1, a typical nu- 
merical result for the real part of the ac conductance for p = 0.52 in a 20x20 size 
sample with leaky capacitors. There may be some genuine concerns whether the 
upper saturation of the ac conductance does actually occur in reahstic samples. In 
response we note that in most of the early experiments, the steady frequency re- 
quired to approach the upper ac saturation regime were probably too high to be 
accessible. Further, we find at least one experiment on Li-doped NiO sample |T^, 
where the upper saturation is clearly observed. 

As a pedagogical example (which will be useful for our interpretations of the 
EM A and the simulation results in the sequel), one may calculate the complex 
G{uj) for some simple prototype circuits using the capacitive conductance gt = jujc 
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where c is the microscopic capacitance and j = a/— 1. For the case of Fig. 2(a): 



On the other hand, for the case of Fig. 2(b): 

It may be noted that the dc conductance G{uj = 0) of the first circuit is l/(ri + r2) > 
(thus it is a conductor), whereas the dc conductance of the second circuit is zero (an 
insulator). The extremely low frequency {uj —>■ 0^) behaviors for both the circuits is 
[Re G{lj) — G{uj = 0)] oc and Im G{uj) oc u. Thus, Re G{lj) or Im G{uj) at very 
low uj cannot distinguish between the two types of circuits. But in the same limit, 
[|G(a;)| — G{uj = 0)] oc cj^ for the elementary percolating circuit of Fig. 2(a), while 
[(^(ci;)! oc uj for the elementary non-percolating circuit of Fig. 2(b). As uj oo, the 
upper ac saturation value for the Fig. 2(a) is 1/ri and that for the Fig. 2(b) it is 
1/(^1 + r2), both of which are finite because the capacitor c does not geometrically 
extend from one electrode to the other. 

Note that on adding together many such elementary circuits (i.e., on increasing 
L), the rational algebraic function type behavior of Re G{uj) obtained from Eq. (|I|) 
and Eq. (|2|) changes over to a sigmoidal type function as shown in the Fig. 1 and 
looks qualitatively very similar to the nonlinear dc conductance as a function of V 

0]. Here, Im G{uj) is zero for both u = and oo, with a broad peak in-between. 
An identical fitting function as the one for G{V) in I, fits (solid line) the simulation 
data (open circles) very well for six decades. In the same spirit, |G(ti;)| in the ac 
case is writen as: 

\Gicu)\ = Gico = 0) + G,(p)[l - exp(-Ao^^)]^ (3) 

To obtain the power-law behavior, one linearizes the exponential function in Eq. (^) 
for small u (much below the upper ac saturation regime) and one gets: [GKcu)! — 
G{0)] oc uj°', where a = fi'y. Now, for considering the ac response in our GRG model. 



we have virtually a three component mixture of the ohmic bonds (conductance 
Qo = g), the t-bonds (in general leaky with a complex gt = Qc + j^c), and the 
insulating bonds {gi = 0). We assume that a sinusoidal voltage V = Voexp{jut) is 
applied across the network. We set Vq = 1, c = 1 and go = ^ for convenience, thereby 
setting the scales for the voltage, the frequency and the conductance respectively. 
We shall concentrate on either the Re G and Im G, or the IGI as a function of u. 
We would also study the phase-angle of the complex G relative to the phase of the 
voltage source at any time t. 

The rest of the paper is organized as follows. In the Sec. 2, we discuss the 
effective medium approximation {EM A) results for our CRC model. Next in the 
Sec. 3, we present the numerical Kirchhof's law solution for the same model and 
compare them with the EM A (for the CRC) and some basic experimental results. 
Finally, in the Sec. 4, we summarize our findings with some concluding remarks. 

2 The EMA result for the CRC Model 

For studying the ac conductance of the correlated RC model within the EMA, we 
assume for simplicity that the t-bonds behave as perfect capacitors [gt = j(^c). So, 
we are in the lower linear dc regime. The probabilities of a bond to be ohmic, 
tunneling (capacitive) or purely insulating in a 2D square lattice are: 

Po = P, (4) 
Pt = ip^ + 3p\ + 3pqYq, (5) 
Pi = l-Po-Pt = [l-ip' + Sp^q + 3pqY]q, (6) 

where q = 1 ~ p. The EMA equation for the effective (i.e., average) complex 
conductance Gg for this three component system may be written as |0, |15[: 

AGl + BGe + C = 0, (7) 
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where A = {d-l^, B = {d-l)[{l-dPo)go+{l-dPt)gt] and C = -[(d-l)-dP,]gogt ■ 
The solution of the quadratic Eq. (|^) is 

-B ± {B' - AACy/' 

Ge = ^ . (8) 

We show below that only the '+' sign in front of the square-root is the physical 
one for our purpose and henceforth use that only. In passing we note that the new 
percolation threshold for our RRTN on a square lattice is Pct = 1/4 using the EM A 
(see I) and pd — 0.181 pl using a real space renormalization |l6l with finite size 
scaling analysis. 

Separating the real and the imaginary parts of Ge{uj), one gets 

Re GM = ^^^^ + ^[X + (X^ + Y^r', (9) 

and, 

Im GeH = ^^^^^f^ + ^[-^ + + ^')'^']'^'' (10) 
where X = {2Po - 1)^ - uj\2Pt - if and Y = 2uj[{2Po - l){2Pt - 1) - 2{2Pi - 1)]. 
It may be noted here that Pt{p) has a single broad peak structure with a maximum 
value of about 0.3840 at a p = 0.4800. Thus, the quantity uj{2Pt — 1) is always 
negative. Further, X is also negative for uj > \{2Po — l)/{2Pt — 1)| and approaches 
— oo quadratically as u; ^ oo. We will take only the absolute value of the square- 
rooted expression for two reasons; (i) Re Ge{uj) in Fig. 3 achieves the necessary 
upper ac saturation since the built-in square-root function in the computer uses 
exactly that, and (ii) this procedure keeps Im Ge{uj) > for all a; > (needed since 
Im gt > 0). 

Let us first check the a; — >■ oo limit. The dissipative part of the complex conduc- 
tance. Re Ge{uj), should be positive, finite, and greater than Ge{uj = 0) in this limit. 
Further, the reactive part Im Ge{uj) of the GRC network must become zero (i.e., 
show no response) when the driving field oscillates much faster than the network's 
relaxation time (or, the time- constant). Now, in this limit. 

Re GM = ^^^^ + I l^^f^ - i^Po - 1)1, (11) 



and, 

Jm G,{u) = ^^i^^LZ}! + 1 |^(2P, _ 1)1, (12) 

and thus for all p > pd, Ge(0) < Re Ge{uj = oo) < oo. Also, clearly Im Ge{uj) in 
this limit is zero. Thus, both the conditions hold in the uj ^ oo limit. 

Next we check the asymptotic expansion as — > 0^. In this limit, [Re Ge{uj) — 
G'e(O)] ~ uj"^ and Im Ge{uj) ~ uj. We find that this extremely low-cu behavior of the 
complex Ge in the EM A is generic for all pd < p < I {except a.t p = pc, see below), 
and [|G'e(^^)| — Ge{uj = 0)] ~ uj"^. Hence the EM A cannot distinguish between a 
percolating and a non-percolating configuration from their low-cj behavior. 

Now, in the special case when p = Pc (= 1/2 for a square lattice), we find that 
in the limit uj —>■ 0^, both the Re Ge{uj) and Im Ge{uj) varies as uj", where a = 
0.5. Obviously, |G'e(t^)| ~ ^^"'^ in this limit. In passing we would like to quote the 
extremely low-frequency EM A exponent in the case of 3D. By looking at the EM A 
expressions for Ge{uj) which are the analogues of Eqs. (^) and (|T0|) for a simple 
cubic lattice at its EM A percolation threshold {pc = 1/3 in 3D), one finds again 
that a = 0.5 (in 3D). 

We show in Fig. 3 a log-log plot of Re Ge{uj) against u in 2D [Eq. (|^)]. In 
conformity with the asymptotic expansions obtained above analytically for very 
small cij's, the EM A results shown in the Fig. 3 give a{p) = 2.0 for all pd < P < 1-0 
except for the special case of p = pc- We have a{pc) = 0.5 both in the very low and in 
the moderately low-u; regimes much below the upper saturation of Re Ge{uj = oo). 
Further for each fixed p ^ pc, there is a characteristic uJo{p) around which a{p) starts 
crossing over from 2.0 to the moderately low-cj exponent of about 0.5. The jump 
of a{p) at an extremely low-cj from 2.0 to 0.5 is the hallmark of the inadequacy of 
EM A. One does also note that for a fixed p Pc , this crossover region becomes 
smaller for smaller p's. Indeed, the crossover region finally tends to vanish as p — > pd 

0.18 in 2D). Thus, we observe that there are no low to moderate u crossovers 
for for p = Pc with a = 0.5 and similarly for p = pd where a = 2.0. 

We did also calculate the EMA results with leaky capacitors at each t-bond 
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{dt = dc+j^^c), i-G., in the upper linear dc regime. The only change here compared to 
the case of perfect capacitors above is that even atp — Pc, |Ge(<^)| — Geij^ — 0)+A;a;^ 
(A;=constant) for very low-o;. Also for all p > Pct, Ge{oj = 0) > 0, as expected. 

We would also hke to see the behavior of phase-angle of the complex conductance 
Ge{uj) with respect to the frequency (a;). The phase-angle (0) is defined through 

ImGeiuj) 

tan(/)=— — — -. 13 
He Ge[u)) 

In the Fig. 4, we plot the phase-angle (0) of the complex conductance against fre- 
quency (uj) for the EM A. As expected, this angle is zero [just like the Im Ge(<^)] 
both at very small and at very large cj's. Further the phase has a peak value = (p^ 
which increases as p is decreased and the lu at which the peak occurs is p-dependent. 
We find from this Fig. 4 that the positions of the peaks tend towards zero and that 
(pm becomes progressively larger as p approaches Pc from higher values. It may be 
noted that we cannot calculate the phase angle for p < Pc in the EM A with 
confidence because the quantity Ge{uj = 0) = 2Po — 1 takes on unphysical negative 
values. Hence we have not shown any curve for p < p^ in Fig. 4. 



3 The Numerical Result for the CRC Model 

We now solve Kirchoff 's laws in our 2D complex network at each node of our corre- 
lated RC model. We obtain the complex conductance of the macroscopic samples, 
their real and imaginary parts, the modulus values and the phase-angle through it- 
erative numerical solution using the Gauss-Seidel relaxation method. In the Fig. 5, 
we have plotted the modulus of the average complex conductance, | < G{uj) > \ 
against u for 0.3 < p < 0.7 and for an external sinusoidal voltage V = cosiujt). 
We let the i-bonds be leaky capacitors with a conductance Qc — 0.001 « Qo and 
averaged over the same twenty configurations for each uj. So, we are virtually in the 
upper dc saturation regime. For clarity, we have shown the graphs in Fig. 5 from 
zero to a moderately low uj (=0.1) much below the upper ac saturation regime. We 
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find that for p , the respective graphs in Fig. 1 and Fig. 5 are fitted very well 
with the Eq. (^) such that a = /X7 ~ 0.7. Thus, in the moderately low-cj regime, a 
obtained experimentally in varieties of systems |Tl], |T^, and in an 'Extended Pair 
Approximation' (EPA) theory is close to what we obtain here (i.e., 0.7). We 
remark here that for a 2D system exactly at Pc , a = 0.5 both in the simple RC 
model (see e.g., Ref. 0) and the EM A of our correlated RC model. Indeed, both of 
them fall short of the realistic value of a not only a.t p = pc, but also at any p ^ Pc- 
Further, these theoretical and experimental results do not support any evidence for 
the statement that a = 6 [|lOl- For example, for our RRTN model near p^ in 2D, 



5^1, whereas a ^ 0.7. Indeed, as we had commented recently [|T^, there is no 
good physical argument supporting the equality. The main point in the comment 
is that since the actual experimental data become flat at large u, it is quite well- 
known that the exponent obtained may depend crucially on the arbitrarily chosen 
(i.e., without any particular physical reason) upper cut-off in u for fitting purposes. 
Consequently, the possibility of an unique function should be explored (as we do 
here) which fits all the way from the lower to the upper saturation range (if suffi- 
cient data are available), and then the power-law behavior in the moderately low-u; 
range should be extracted from that function. 

Now, as we had mentioned before, we know of at least one experiment by Pollak 



et. a/.|[Tj] on Li-doped NiO single crystals at sufficiently low temperatures, where 
the upper saturation of the Re Ge{uj) may be clearly observed. The frequency range 
used in that work is from about 10^ Hz to about 10^° Hz. We did not try to fit 
them by our method since no data below 10^ Hz was available. In any case we 
note that the upper saturation is also consistent with the fact that the measured 
relaxation time for this sample is about 2.2 x 10~^° s. In passing we do also note 
that in this system as well, the doping concentration of Li is extremely low: from 
about 13 X 10~^ to about 136 x 10^®, and that the pristine NiO is an insulator. We 
observe that the upper saturation value of the conductivity for this experiment [|1^ 
is about 0.1 S/m, and that the data for all the three concentrations used for Li, fall 
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closely enough to their fitting function in the high frequency range, but not so well 
in the low/ moderate frequency range. The power-law exponent a (around VS^Hz) 
seems to be close to 2.0. 

Next, we we would like to emphasize that the apparently excellent fitting shown 
in Fig. 1 with Eq. may be misleading to the eye in the very low-cj regime. A 
careful analysis of data in a more revealing log-log plot shows that for very low cj, 
the fittings may be actually quite bad. We observed for the two elementary circuits 
(in Fig. 2) for very low cj, |G'(c(j)| ~ cj^ or depending upon whether the circuit is 
percolating or not. Similarly, we anticipate that for extremely low lo and in the lower 
dc regime, the CRC network would also show such simple behaviors (instead of a 
more complicated percolative behavior in the moderately low-cu regime as described 
above). Further this simple behavior is expected to persist (as seen in our EM A 
results) generically for each p up to some scaled crossover frequency ojq depending 
on Qo, Qc and c. For oj > Uq, we expect macroscopic percolative effects to gain 
control and instead of Eq. (H), |G'(co')| should follow a more general equation closer 
in form to that used for the dc- conductance [cf. Eq. (4.7) of I]: 

1^(0;) I = G{uo) + G,{p)[l - exp(-A[^ - uof)?- (14) 

A typical fit by the above equation is shown in Fig. 6 for a system size L = 20 
and p = 0.3 for about seven decades in u using leaky capacitors. In the inset of 
Fig. 6, we show the log-log plot for the modulus value of G and observe a very 
convincing quadratic behavior with | < G{uj) > \ = 0.01219 + 55u;^ for ujq < 0.01. 
This observation is totally matching with our analyses for Fig. 2(a). For u > uq, we 
get a very good fit for the intermediate frequency range with the general Eq. (|T^) 
where a = /i7 ~ 3.0 (which is much larger than 0.7). 

Next in the Fig. 7, we show another fit for a single configuration (at p = 0.45 < 
Pc) with perfect capacitors at the t-bonds {qc = 0, non-percolating) of the CRC 
network with L=20. In this case one can easily observe that in the very low-cj range 
up to a crossover frequency uq — 0.01, |G'(co')| ~ Gu. This behavior is akin to that 
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of the elementary circuit of Fig. 2(b). Further, beyond ujq , we do again have an 
excellent fitting with the Eq. ([T^) with a = fi'y c::^ 0.5. So we have two things to note 
here for non-percolating configurations (with perfect capacitors, i.e., in the lower 
dc regime): (i) the very low-u behavior is linear in uj, and (ii) the intermediate 
frequency behavior seems to give a lower value for the exponent compared to that 
for percolating configurations (e.g., as in the upper dc regimes of Fig. 1 and Fig. 6 
where 6 = 0.7 near pc). Note that for the non-percolating configurations (only at 
Pc) in the lower linear dc regime, a = 0.5 in the case of EM A as well. 

For many practical situations, this intermediate frequency range (away from both 
the lower and the upper ac saturation regimes) is of prime interest. For t-bonds with 
leaky capacitors at any p > pd (i.e., in the upper dc regime), we find that a has a 
minimum value of about 0.7 near pc, and increases on both sides of it with a value of 
about 3.0 at p = 0.3 (as shown above) and of about 1.35 (not explicitly shown here) 
a.t p = 0.7. Clearly, this result has a qualitative similarity with the dc nonlinearity 
exponent 6{p) as a function of p as shown in the Fig. 11 of I. 

Finally, the variation of the phase-angle (0) with frequency {u) obtained by 
numerical (Kirchoff's laws) solution of our CRC network has been shown in the 
Fig. 8 for 0.3 < p < 0.7. One can easily observe from this figure that in the CRC 
model, for configurations with p around Pc for a 2D square lattice, the peak value 
of the phase, 0^ — 0.7 (radian). This is close to the universal phase-angle value 
of 7r/4 radian obtained in the simple RC model in 2D at pc as predicted by Clerc 
et. a/. 0. As noted before, the EM A calculations (i.e., mean- field approach) of the 
phase (piuj) for our CRC model also obtains 0m = 0.7 rad when p = pc- The main 
problem with the EMA in this case is that the peak occurs at an u of about 0.1 
rad/s (instead of at around 0.2 rad/s as for the CRC in Fig. 8 above) and that it is 
too broad. Further, as shown in the Fig. 8 for the CRC network which is expected 
to be much more realistic (better than the mean-field theory, for one reason) from 
the demonstrations and the arguments above, 0m increases as p decreases: from a 
value of 0m — 1-1 rad for p = 0.3 to a value of 0m — 0.4 rad for p = 0.7. 
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4 Summary and Conclusion 



In this work we have studied the nonhnear ac response (in the two hnear dc regimes) 
of a variety of materials, e.g., a binary metal-insulator type composite, a dispersed 
metal, amorphous semiconductors at high field, etc. by modeling them in terms of 
our recently proposed RRTN |^, |^ model. This RRTN (Random Resistor cum 
Tunneling-bond Network) driven by an ac electric field gives rise to a correlated RC 
{CRC) model because we place the capacitors only at the t-bonds whose positions 
are correlated to the positions of the nearest neighbor gaps between two ohmic bonds, 
as discussed above. For simplicity, the rest of the insulating bonds are assumed to 
remain inactive at any voltage however large. Thus this model effectively behaves 
like a three component mixture. We studied the complex ac conductance G{uj) 
of our CRC model as a function of uj and at a fixed low voltage both by using an 
effective medium approximation {EM A) and by an iterative numerical solution using 
Kirchoff's laws. The real part Re G{uj) or the modulus |G'(co')| behaves qualitatively 
very similarly to that for the real, nonlinear dc conductance ^(V") as a function 
of the external voltage (see I). Thus we fit the Re G{uj) or the |G'(t^)| with an 
exponential- related function as in Eg. ([T^) . The fitting with this function seems to 
be very good just like the fitting of G{V) with a similar function is optimum. But 
the similarity ends there. 

Whereas there is only one exponent for the power-law type low voltage dc re- 
sponse, there are two different power-law regimes for, say [^^(u;)!, much below its ac 
upper saturation regime. Since any realistic sample should not be able to respond to 
an extremely large uj, Im G{u!) must vanish a.s uj —>■ oo, and that does automatically 
happen in our CRC model. Thus \G{uj = oo) \ is real and finite. For non-percolating 
systems at very low-tu, |G(co')| ~ u, whereas if percolating, [[^(co')! — G(0)] ~ cu^. 
This happens irrespective of the value of p. Next for uj > ujq, i.e., at moderately low 
tu's (still much below the upper ac saturation regime), the system crosses over to 
a non-integer power law behavior with a positive exponent a{p) which depends on 
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the value of p. If we use leaky capacitors at the t-bonds, the system would percolate 
for all p > Pet (— 0.18) 0. For such situations, if p ~ pc, a has the minimum value 
of about 0.7 and this value matches with that for many experiments. Away from 
Pc (both for higher and lower values), with leaky capacitors at the t-bonds, a{p) 
becomes progressively larger than 0.7. Now, if we use perfect capacitors, then for 
non-percolating configurations close to pc, <y{p) — 0.5 which is significantly lower 
than 0.7. It may be noted that both the EM A and the traditional RC model cannot 
obtain anything other than the value of 0.5 for p close to pc- Further, in the case of 
EM A, the exponent a{p) = 2.0 for any p ^ Pc- 

We did also look at the phase angle 4>{uj) of the complex G at several p's. Since 
Im G{uj) vanishes both at very high and very low cu's (at any p), (j){uj) would also 
have this property as a function of u with at least one peak value (0m) somewhere 
in-between. There is one curve for each p with a different peak value. In the case 
of EM A, (f){uj) cannot be calculated for any p < Pc and so in Fig. 4, we have shown 
them only for p > p^. For the CRC network, we do not have this pathology, and 
hence in the Fig. 8, we have shown it for 0.3 < p < 0.7. In general, 0^ increases 
both in the EM A and our numerical solution of the CRC model as p decreases. For 
the latter, it is about 0.4 rad aX p = 0.7 and 1.1 rad for p = 0.3. It is interesting 
to note that at p = Pc, (pm takes the universal value of about 7r/4 rad in the simple 
RC model, numerical solution and the EM A of our CRC model. But, whereas the 
4'm{Pc) occurs at an C(j = 0.2 rad/s for the CRC model, it occurs at an u; = 0.1 rad/s 
in the case of the EMA. Further the peak around 0^ is too much broad for the 
EMA. 
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Figure Captions: 

Fig. 1 Real part of the complex G{uj) for our model correlated RC (CRC) 
network a.t p = 0.52, with leaky capacitors placed at the tunnel junctions. The solid 
line given by Eq. (j^) appears to give a very good fit. 

Fig. 2 Two prototypical elementary circuits each with two ohmic resistors and 
one capacitor. Note that the Fig. 2(a) corresponds to a metal-like zero-frequency be- 
havior (i.e., Ge{uj = 0) > 0), whereas the elementary circuit in Fig. 2(b) corresponds 
to an insulator (i.e., Ge{uj = 0) = 0). 

Fig. 3 The EM A result in 2D for the CRC model: real part of the effective 
conductance Re Ge{uj) against uj for a set of values of p. 

Fig. 4 The EM A result for CRC model in 2D: the phase-angle (0) against uj 
for different values of p. 

Fig .5 The numerical results in 2D for the CRC model with leaky capacitors 
with a finite real conductance Qc > 0: the | < G{uj) > \ against uj for a set of values 
of p. For each p and u, the average over 20 different configurations were taken. 

Fig. 6 The | < G{uj) > \ for p = 0.3 and a square lattice of size L = 20. 
Just as in Fig. 5, each t-bond represented a leaky capacitor and an average over 
20 configurations were taken. The inset shows an tu^-dependence upto a crossover 
frequency of ujq ~ 0.01. Above ujq, the Eq. (jl^) gives a very good fit as shown. 

Fig. 7 Another example of |G'(c(j)| against u for p = 0.45 and for a typical 
configuration on a square lattice of size L = 20. The very low-cu part shows a purely 
linear behavior. But for a; > cuq — 0.01, Eq. (|I^) gives the optimum fit. 

Fig. 8 Phase-angle (0) as a function of u from the numerical solutions of our 
CRC model. 
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